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In this article, the reliabilities R(t) = P(X > t), when X follows two-parameter geo¬ 
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1 Introduction 


Various lifetime models have been proposed to represent lifetime data. Most of these models 
assume lifetime to be a continuous random variable. However, it is sometimes impossible or 
inconvenient to measure the life length of a device on a continuous scale. In practice, we come 
across situations where lifetimes are recorded on a discrete scale. Discrete life distributions 
have been mentioned by Barlow and Proschan pQ. Here one may consider lifetime to be the 
number of successful cycles or operations of a device before failure. For example, the bulb in 
xerox machine lights up each time a copy is taken. A spring may breakdown after completing 
a certain number of cycles of ‘to-and-fro’ movements. 

The study of discrete distributions in lifetime models is not very old. Yakub and Khan [2] 
considered the geometric distribution as a failure law in life testing and obtained various para¬ 
metric and nonparametric estimation procedures for reliability characteristics. Bhattacharya 
and Kumar [3] discussed the parametric as well as Bayesian approach to the estimation of the 
mean life cycle and reliability for this model for complete as well as censored sample. Krishna 
and Jain [4] obtained classical and Bayes estimation of reliability for some basic system con¬ 
figurations. Modeling in terms of two-paranreter geometric and estimation of its parameters 
and related functions are of special interest to a manufacturer who wishes to offer a minimum 
warranty life cycle of the items produced. 

The two-paranreter geometric distribution [abbreviated as Geo(r, 6)] given by 

P(X = x) = (1 -d)9 x ~ r ; x = r,r + l,r + 2,... 0 < 9 < 1, r € {0,1, 2,...}, (1.1) 

is the discrete analog of two-parameter exponential distribution. If X follows a two-parameter 
exponential distribution, [Y], the integer part of X, has a two-parameter geometric distribution 
[see Kalbfleish and Prentice [5l Ch. 3]]. The reliability of a component when X follows two- 
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parameter geometric distribution is given by 


R(t) = O t ~ r ; t = r,r + l,r + 2,.... (1.2) 

Laurent [U] and Tate [7] obtained the uniformly minimum variance unbiased estimator (UMVUE) 
of the reliability function for the two-parameter exponential model. Different estimators of this 
reliability function have been discussed in Sinha [8], 

If a system consists of m identical components each follows two-paranreter geometric distribu¬ 
tion, then the reliability of fc-out-of-m system is given by 

R s (t) = P(X {m _ k+1) >i) = J2(™)R(tY[ l-m)}™-*; t = r,r + 1, r + 2,... (1.3) 

Special cases of R s (t) give series system (for k = m) and parallel system (for k = 1). 

In the stress-strength setup, R = P(X < Y ) originated in the context of the reliability of 
a component of strength Y subjected to a stress X. The component fails if at any time the 
applied stress is greater than its strength and there is no failure when X < Y. Thus R is a 
measure of the reliability of the component. Many authors considered the problem of estimation 
of R in continuous setup in the past. Particularly, for the two-parameter exponential set up, 
Beg m derived the MLE and the UMVUE of R. In the discrete setup, the reference list is 
very limited. Maiti m has considered stress (or demand) X and strength (or supply) Y as 
independently distributed geometric random variables, whereas Sathe and Dixit m assumed 
as negative binomial variables, and derived both MLE and UMVUE of R. Maiti and Kanji 
m has derived some expressions of R using a characterization and Maiti [T3TT41 considered 
MLE, UMVUE and Bayes Estimation of R for some discrete distributions useful in life testing. 
All the above mentioned works have been concentrated on one-parameter family of discrete 
distributions. 
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If X and Y follow two-parameter geometric distributions with parameters (■ Oi,ri ) and 
(O 2 , r-i) respectively, then 

R = p9 2 for S > 0 

= 1-(1 -p)0i S for d < 0, (1.4) 

where p = and S = r\ — ■ 

Here we are interested to see whether similar estimators are obtained in case of the two- 
parameter geometric distribution, the discrete analog of the two-parameter exponential dis¬ 
tribution. Then, it might be straightforward to use the two-parameter geometric distribution 
in the discrete life testing problem where a minimum warranty life cycle of the item is offered. 
In this article, we have found out some estimators of both R(t) and R s (t) for complete as well 
as censored sample. Some estimators of R have also been provided. The estimators have been 
compared through simulation study. 

The paper is organized as follows. In section 2, we have derived MLE and UE of both R(t) 
and R s (t). We have also derived MLE of these reliability functions for type-I censored sample. 
MLE and an unbiased estimator of R have been found out in section 3. In section 4, simulation 
results have been reported. Section 5 concludes. 

2 Estimation of R(t) and R s (t) 

Let (Xi,X 2 , ...,X n ) be a random sample from Geo(r,9) and (XL), X/ 2 ), ...,X( n )) be ordered 
sample. Maximum Likelihood Estimator of r and 9 are X^ and respectively, where 
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S = Y^i =i (Xi — -X’(i)) = Y^=i (X(i) — -X’(i)) . ML Estimators of R(t) and R s (t) are given by 


R M (t) = 1 for t < 


S 


n + S 


t-x, 


(i) 


for t > X, 


(i) 


(2.5) 


and 


RsAi(t) = 1 for t < X( X ) 


i=k 


^2 ( • ) RM(t) l - Rm(R) 


for t > X, 


(i) 


( 2 . 6 ) 


respectively. 

Suppose we record the observations (X( 1 ),X( 2 ), ...,X^),p < n that are failed before a pre¬ 
specified number of cycles c and remainings survive beyond c. Then, MLE of r and 9 are X^) 
and respectively, where S* = X)?=i ~ ^(i)) + (n — p){(c + 1) — X^j}. Hence ML 

Estimators of R(t) and R s (t ) are given by 


R*M{t) = l 


for t < X(x) 


S* 


p + S* 


t-x , 


(i) 


for t > X, 


(i) 


(2.7) 


and 


R*sM{t ) = 1 


for t < X(x) 


i=k 




for t > X, 


(i) 


( 2 . 8 ) 


respectively. 


Theorem 2.1 (X^S) is sufficient statistic for ( r,9 ). 


Proof: Let X = (Xi, X 2 ,..., X n ), we have to prove that P (X = x|Xo) = u, S = s) does not 
depend on r and 9. 


5 



Given = u and S = s, X_ is an n-dimensional random variable with domain A UjS = 

{x\x {1) = u, Ya= 2 x i = s + (n - 1 )u}. 

For x € A U:S , 

and 

n 

p(x = x) = JJ p(Xi = Xi ) = (i - e) n e s+n{ ~ u ~ r \ 

i= 1 

Thus P(X = x\X^ = u, S = s) = i ^ 1 | , where |^4 U)S | is the number of elements in ^4 UjS . The 
number of elements in A UjS is the number of possible n-uplets (x\, x%, ..., x n ) such that = u 
and Yl'i =2 x i = s + (n — 1 )u which clearly does not depend on 9 and r. 

But (A - ^, S) is not complete as it is to be seen from the following counter example. 

Counterexample 2.1 Let us define g(., .) as 

g{X (1) ,S) = 1 if X (1) = r + 2, 5 = 0 

= —1 ifX(i)=r + 1, S = n 
= 0 otherwise. 

Now A'(x) and S can take values r + 2 and 0 with the probability (1 — 9) n 9 2n (for X( 2 ) = 
r + 2 , Xr n \ = r + 2), and A^) and S can take values r + 1 and n with the probability 
(1 — 9) n 9 2n (one such particular situation is Xr 2 ) = r + 2 , Ap^) = r + 2 , Xr n \ = r + 3). 
Therefore, it is found that E T) g [g (Xm, 5)] = 0 but g (Am, S ) 7 ^ 0. 

The upcoming theorem will demonstrate the conditional distribution of X for given (Am, S'). 

Theorem 2.2 The conditional distribution of X given (A^pS) is as following: 

For n = 1, 

P{X = x/X w ,S) =1 


6 



For n = 2 , 


P(X = x/X (1] ,S) = 1 if x = X (1) 

= \ if x = X ( i) + 5 

For n > 3, S < n, 

r(x =x/ x m ,s) - >fx m <*<x m + s 

= 0 otherwise. 

For n > 3, S > n, 

. rr ! )/{( s *r')-C:0) *■-*<» 

ifX (1) <x<X (1) +S-(n- 1) 
/5-(x-X ( 1 ) )+n-2\ f/S + n-l\ _ (S - A 1 

V s-(*-x (1) ) H ^ y U-vi 

if + 5 — (n — l)<x< X(!) + S 

= 0 otherwise. 

Proof: Joint distribution of X \, X 2 , X n is given by 

P(X 1 , X 2 , X n ) = (1 - 6») n 0^=i( Xi - r ) = (1 - 0 )" 0 S+n(* ( 1 )-r) } r < X (1) . 

Now, 

p , y _ _ Ex,, x, .X, PjX = *. V2. .... X n /x m , s ) 

1 j/ (i), ) Exi, x„ x„ .... x„ P(Xu x 2 , XJX m , S) 

E(X,, X, . x„/x (1) , s) 1 

E(X1, X,, X,. X„/X (1) , 5) 1 

Here the denominator is equivalent to finding out the total number of ways in which S indis¬ 
tinguishable balls can be placed in n cells so that at least one cell remains empty. In general, if 
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there are r indistinguishable balls to be placed randomly in k cells, then the number of distin¬ 
guishable distributions is whereas the number of distinguishable distributions in which 


no cells remains empty is Therefore, we get the denominator as ( S+ g X ) — (f_i) if S > 


n 


and if S < n, it will be ( + 5 _1 )- Similarly, the numerator is equivalent to finding out the total 

number of ways in which S — (x — Xm) indistinguishable balls can be placed in n — 1 cells so 

that at least one cell remains empty and hence, we get it as ( 5 £ r 2 ) — ( S *) 

v o — {%— A^(i) ) 7 v n—z / 

if S — (x — X(i)) > n — 1 and if S — (x — < ra — 1, it will be ( 5 ^ 2 )' Hence the 

theorem follows. 

Since S) is sufficient but not complete statistic for (r, 6), we are handicapped of searching 

the UMVUE of any estimable function of these parameters using the Lehmann-Scheffe theorem. 
Hence, we will find an improved estimator of the reliability functions using the Rao-Blackwell 
theorem. 


Define 


Y = 1 if X, > t 

= 0 otherwise. 

Then R{t) = E(Y ) = P (X\ > t). Using the Rao-Blackwell theorem, an unbiased estimator of 
R(t) is given as follows: 
for n = 1 , 

Ru{t) = 1 if t < 

= 0 if t > X(iy 
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for n = 2 


Ru(t) = 1 if t < 


if X 


(1) < t S A (1) 


+ 5 


= 0 if f > X (1) + 5. 


for n > 3 and 5 < n, 
Ru(t) 


= 0 


1 

if t < x (1) 


A (1 )+S 

v 

/5 — (z — 

W(i)) + n — 2 

x=t 

V 5-1 

! x - x <i)) ^ 

0 

if t > X (1) + 5. 


S 


if X (1) < t < X (1) + 5 


It can also be written as 


Ruit) = 


1 

A ( i)+5 

E 


X=t 


= 0 


if t < x (1) 

n — 1 

X (i) + S' + n - - j=1 
if t > X (1) + 


n —2 


in 


(X(i) + S + n-x-l-j) . 
(X(i) + S + n — 1 — j) 


if X, 


(1) < t S a (1) 


+ 5 


For n > 3, S > n, 
Ruit) = 1 


if t<X, 


^ A(i) 


X m +S-(n- 1) 


E 


S-(x- X (1) ) + n - 2\ _ fS - (x - X (1) ) - 1 


*=* l V s-(x- x (1) ) 

A(i)+5 

x =X (1) +S—(n—1)+1 ^ S - [x- X(!)) 


n — 2 



5 + n — 1 
5 


5 - (z - X (1) ) + n — 2\ J/5 + n — 1 



5 


5-1 
n — 1 


if X (1) < t < X (1) + 5 - (n - 1) 


- (z - X (1) ) + n-2\ ((S + n-\ 

h v ^(z-x^ ]/ 



s 


5-1 
n — 1 


if X (1) + 5 - (n - 1) < t < X (1) + 5 


= 0 


otherwise. 
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In other way the estimator Rjj(t) is to be UMVUE if it is uncorrelated with all unbiased estima¬ 
tor of zero. We take a class of unbiased estimator of zero as Uq = {u : y” =1 qA* = u, Yli =l c * = 
1}. If Rjj(t ) is UMVUE, then Cov(Uq, Ru(t )) = 0 i.e. Ccw(1000.t/o, WOO.Ru(t)) = 0. Ana¬ 
lytical derivation seems to be intractable. We go for simulation study taking some particular 
choices of (ci, C 2 , ..., c n ) and different t, and 1000 covariances have been calculated and their 
averages have been shown in Tables 7-8. It is noticed that they are not uncorrelated and hence 
Ruit) is not UMVUE. 

y I(X>t) 

The variance of this unbiased estimator will be smaller than the unbiased estimator — !=1 — . 

n ’ 

where /(.) is the indicator function. 

To study the asymptotic behavior of Rjj(t) we conduct a simulation study taking different 
values of parameters. 10000 estimates of Rjj(t ) and R.M{t), their variances, 95% confidence 
limits and coverage probability (CP) have been shown in table 9. Histogram of Ruit ) for 
n = 20, r = 15, t = 25, 8 = 0.96 has been shown in Figure |T| In this set up the true reliability, 
R(t) = 0.6648326. The figure is near normal. From the table 9, it is also evident from coverage 
probability point of view, Rm( t) is better if 0.02 < R(t) < 0.5, otherwise Ruit ) is better. From 
the table it is observed that asymptotic variance is approximately _ 

Define 

Z = 1 if at least k of X^s among X±, X 2 , ..., X m are greater than or equal to t 

= 0 otherwise. 
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Figure 1: Histogram of Ruit). 

Then R s (t) = E(Z) = (7)[-^(0P[l ~ R{t)] m ~ l . Using the Rao-Blackwell theorem, 

unbiased estimator of R s {t) is given by (for 2 < m < n) 


R s u{t ) 


1 if t < X(!) 

£ (j) [A/Wl‘11 - Ru(t)} m ~‘ if X m < t < x m + s 

0 if t>X (1) + S. 


an 
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3 Estimation of R 


Let (Xi, X 2 ,X ni ) and (Yi,Y 2 , ...,Y n2 ) be random samples from Geo(r 1 ,0 1 ) and Geo(r 2 ,0 2 ) 
respectively. (XWjSi) and are defined in the same way as in section 2. Hence ML 

Estimator of R is given by 

for 5 > 0 

= 1 - (1 - P) ( x Q ) for <5 < 0, 

\m + Si) 

where p = and 8 = X {1) - Y (1) . 

We define censored scheme in the same way as in section 2, with pre-specified censored 
cycles ci and C 2 and with p\ and P 2 censored observations. 

Then, ML Estimator of R is given by 

”’ M = '’’Gifs-) fo, '’ >0 

= 1 “ (1 “ h Grfsr) for,s<0 ’ 

where P* = and * = X (D ~ Y W- 

Define 


Rm = P 


S 2 

n 2 + S -2 


Z = 1 

= 0 


if X, < Y 1 
otherwise. 
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Then R = E(Z ) = P (X\ < Yi). Application of the Rao-Blackwell theorem gives an unbiased 
estimator of R as 

Y(i) min(W-L,W 2 ) W2 

Ru = — + E f(x/X ( i),S 1 ) + E E f(x/X (1] ,S 1 )f(y/Y (1) ,S 2 ) if X (1) < y (1) 

ni ®=X (1) +1 x=Y ( i) 

min(W\,W2) W2 

= -+ E E/(*/^ ( i)^i)/(y/y(i)^2) if at ( i) = V(i) 

®=X (1) J/=X 

W2 min(Wi,W 2 ) W2 

= - E f(y/ Y W’ S *) + E '£f(x/X w ,S 1 )f(y/Y (1) ,S 2 ) ifX (1) >y ( i) 

ni 2/=X (1) x=X (1 )+l ?;=£ 

where, 1-Ti = Xaj +Si, W 2 = Y A (i) + S 2 ■ The variance of this unbiased estimator will be smaller 
than the unbiased estimator X^j= 1 1(Xi < Yj). 

4 Simulation Study 

4.1 Discussion on simulation results relating to R(t) and R s (t). 

In order to have an idea about the selection of an estimator between Maximum Likelihood 
Estimator (MLE) and Unbiased Estimator (UE), Mean Squared Errors (MSEs) and hence 
percent relative efficiency using these MSEs have been calculated for R(t) and R s (t). We 
generate sample of size n and on the basis of this sample, calculate MLE and UE. MLEs have 
been calculated for complete as well as censored (type-I defined in earlier section) samples. 
10000 such estimates have been calculated and results, on the basis of these estimates have 
been reported in Tables 1 — 6. Initial set up for parameters has been chosen as (n = 20, r = 
15, c = 25, 9 = 0.8, t = 25, k = 2, m = 8). Each table has been prepared considering different 
choices of a particular parameter, keeping others fixed at initial set up. All simulations and 
calculations have been done using R-Software and algorithms used can be obtained by contacting 
the corresponding author. Different columns of a table are as follows: 
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Col.l: Component Reliability, R(t) 

Col.2: Average of MLEs of R(t) for complete sample 

Col.3: Average of MLEs of R(t) for censored sample 

Col.4: Percent relative efficiency of MLE of R(t) for complete sample to that of censored 
sample 

Col.5: Average of UEs of R(t) 

Col.6: Percent relative efficiency of UE to MLE of R(t) 

Col.7: System Reliability, R s (t.) 

Col.8: Average of MLEs of R s (t) for complete sample 

Col.9: Average of MLEs of R s (t) for censored sample 

Col. 10: Percent relative efficiency of MLE of R s (t) for complete sample to that of censored 
sample 

Col. 11: Average of UEs of R s (t) 

Col. 12: Percent relative efficiency of UE to MLE of R s (t) 

For estimation of component as well as system reliability, MLE performs quite well from 
efficiency point of view except for very less reliable component whose importance in practice 
is not so much meaningful (Table 1 and Table 2, Col.6 and Col. 12). Even though the MLE is 
biased, the combined effect of variance and bias is less than the variance of the proposed UE. If 
sample size n increases, as expected, MLE and UE for component as well as system reliability 
become equally efficient (Table 4, Col.6 and Col. 12). If censored number of cycle c increases, 
MLEs for complete as well as censored sample become equally efficient (Table 3, Col.4 and 
Col.10). 
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4.2 Discussion on simulation results relating to R. 


In order to have an idea about the nature of the MLEs under complete sample and censored sam¬ 
ple, we have calculated percent relative efficiency of Rm with respect to R*m i-e. pf ^ X100 

and presented in tables for n = 5, 10, r2 = 10, 5, 9\ = 0.7, 0.8, 6-2 = 0.7, 0.8 and for different 
values of c\ and C2 (Table 10-13). Here, we take n\ = n-2 = 10, and 1000 estimates of R have 
been taken for calculating MSEs. We observe that, as expected, MLE of R for censored sample 
approaches to that of complete sample as both c\ and C2 increase. Fixing any one of c\ and C2, 
and increasing the remaining, efficiency of MLE for censored sample sometimes increases but 
there is no guarantee. 

In order to have an idea about the selection of an estimator between UE and MLE, we have 
calculated MSEs of Estimates of R. We prepared tables (Tablel4-18) for (i) 9\ = 0.1, 62 = 0.1, 
(ii) 61 = 0.5, 62 = 0.5, (iii) 6 X = 0.8, 9 2 = 0.2, (iv) 6 >i = 0.9, 6 2 = 0.9, (v) 9 1 = 0.2, 9 2 = 0.8,and 
r\ and r2 equal to 5, 10, 15 and 20, n\ = ri2 = 10. In these tables, values in 1st row indicates 
true value of R, 2nd and 3rd rows indicate average of estimate of MLE and UE of R, and 4th 
and 5th rows indicate MSEs of MLE and UE of R respectively. 

We observe that in almost all cases, MLE is better that UE for R in mean square error 
sense. Therefore, as soon as we entered to unbiased class, we are losing some efficiency. It 
is to be noticed that MLE in this case is not an unbiased estimator. Moreover, MLE has a 
computational ease. 

5 Concluding Remark 

This paper takes into account the inferential aspects of reliability with the two-parameter 
geometric lifetime. The continuous distributions are widely referenced probability laws used 
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in reliability and life testing for continuous data. When the lives of some equipments and 


components are being measured by the number of completed cycles of operations or strokes, or 
in case of periodic monitoring of continuous data, the discrete distribution is a natural choice. 
At the same time, if a minimum warranty life cycle of the items are provided, the two-parameter 
geometric distribution is the simplest but an important choice. Under this set up estimators 
of reliability functions- under mission time as well as under stress-strength set up, have been 
viewed. It is interesting to note that, unlike the case of the two-parameter exponential, the 
estimators of the parameters of the two-paranreter geometric distribution are not complete. 
As a result, we only get unbiased estimators of the reliability functions for the two-parameter 
geometric set up. 

In most of the situations, MLE gives better result than the UE in mean square sense. As soon 
as we entered to unbiased class, we are loosing some efficiency. If one is ready to sacrifice the 
unbiased criteria of the estimator, the MLE in this case is preferable. It is to be noticed here 
that MLE is not an unbiased estimator. Moreover, MLE has a computational ease. 
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Table 1: Calculations relating to R(t) and R s (t) 
n = 20, r = 15, c = 25, 9 = 0.8, k = 2, m = 8 



Col.l 

Col.2 

Col. 3 

Cold 

Col. 5 

Col.6 

Col.7 

Col.8 

Col.9 

Col.10 

Col.ll 

Col.12 

16 

0.8 

0.79393 

0.79491 

104.89 

0.80022 

113.64 

0.99992 

0.999708 

0.999706 

100.15 

0.99977 

138.10 

17 

0.64 

0.62917 

0.63338 

106.57 

0.63867 

119.58 

0.99571 

0.982063 

0.982062 

100.31 

0.98545 

128.81 

18 

0.512 

0.49571 

0.50086 

106.73 

0.51097 

109.54 

0.96979 

0.928259 

0.929258 

99.76 

0.93871 

120.28 

19 

0.4096 

0.39881 

0.40459 

110.55 

0.41452 

98.88 

0.90330 

0.836803 

0.838851 

100.28 

0.85352 

111.01 

20 

0.32768 

0.31945 

0.32699 

113.34 

0.33407 

92.92 

0.79549 

0.722729 

0.729116 

101.38 

0.74399 

102.54 

25 

0.10737 

0.10622 

0.11228 

131.55 

0.10647 

80.50 

0.20909 

0.220599 

0.236196 

125.40 

0.22473 

83.46 

30 

0.03518 

0.04018 

0.04422 

151.20 

0.03661 

83.31 

0.03009 

0.056895 

0.067862 

168.06 

0.05466 

83.63 

31 

0.02814 

0.03461 

0.03810 

151.25 

0.03011 

89.96 

0.01981 

0.044496 

0.054024 

178.83 

0.04043 

91.94 

35 

0.01153 

0.01476 

0.01754 

212.54 

0.01108 

104.80 

0.00355 

0.011759 

0.017823 

334.81 

0.00923 

101.20 

40 

0.00378 

0.00692 

0.00868 

210.03 

0.00473 

122.46 

0.00039 

0.003788 

0.006659 

375.59 

0.00271 

107.04 

45 

0.00123 

0.00271 

0.00342 

226.31 

0.00157 

174.54 

4.26 x 10" 5 

0.000858 

0.001667 

737.01 

0.00048 

189.46 































































































































































Table 2: Calculations relating to R(t) and R s (t) 
n = 20, r = 15, c = 25, 9 = 0.8, t = 25, m = 8 


Col.l 

Col. 2 

Col. 3 

Col.4 

Col. 5 

Col.6 

Col. 7 

Col.8 

Col.9 

Col. 10 

Col.11 

Col.12 

0.10737 

0.10841 

0.11121 

116.64 

0.11827 

78.010 

0.59695 

0.56534 

0.57109 

105.37 

0.59365 

94.259 

0.10737 

0.10821 

0.11118 

118.38 

0.11805 

78.110 

0.45796 

0.06392 

0.06981 

137.86 

0.07894 

63.170 

0.10737 

0.10829 

0.11124 

119.96 

0.11813 

78.048 

3.53 x 10" 05 

0.000263 

0.000408 

374.98 

0.000436 

39.115 

0.10737 

0.10933 

0.11229 

118.16 

0.11929 

77.648 

1.76 x 10 -08 

7.12 x 10 -07 

1.51 x 10 -06 

1081.50 

1.46 x 10“ 06 

25.299 


Table 3: Calculations relating to R(t) and R s (t) 
n = 20, r = 15, 9 = 0.8, t = 25, k = 2, m = 8 


c 

Col.l 

Col.2 

Col.3 

Col.4 

Col. 5 

Col.6 

Col.7 

Col.8 

Col.9 

Col. 10 

Col. 11 

Col. 12 

20 

0.10737 

0.10790 

0.11599 

164.86 

0.11772 

78.104 

0.20909 

0.21959 

0.24290 

160.05 

0.24860 

76.145 

25 

0.10737 

0.10819 

0.11123 

118.83 

0.11801 

78.133 

0.20909 

0.22043 

0.22924 

118.50 

0.24941 

76.362 

30 

0.10737 

0.10829 

0.10946 

107.18 

0.11814 

78.029 

0.20909 

0.22059 

0.22397 

107.18 

0.24964 

76.268 

35 

0.10737 

0.10899 

0.10962 

103.95 

0.11892 

77.747 

0.20909 

0.22251 

0.22432 

103.93 

0.25179 

76.101 

40 

0.10737 

0.10842 

0.10866 

101.60 

0.11828 

78.084 

0.20909 

0.22093 

0.22160 

101.54 

0.24995 

76.446 

45 

0.10737 

0.10772 

0.10780 

100.63 

0.11751 

78.329 

0.20909 

0.21914 

0.21935 

100.58 

0.24801 

76.463 
































































































Table 4: Calculations relating to R(t) and R s (t) 
r = 15, c = 25, 9 = 0.8, t = 25, k = 2, m = 8 


n 

Col.l 

Col. 2 

Col. 3 

Col. 4 

Col.5 

Col. 6 

Col. 7 

Col.8 

Col.9 

Col. 10 

Col. 11 

Col. 12 

10 

0.10737 

0.10599 

0.11160 

127.07 

0.12610 

62.804 

0.20909 

0.21957 

0.23384 

122.73 

0.27625 

63.291 

15 

0.10737 

0.10769 

0.11140 

120.18 

0.12087 

72.608 

0.20909 

0.22115 

0.23150 

118.81 

0.25938 

71.594 

20 

0.10737 

0.10909 

0.11298 

118.70 

0.11984 

77.429 

0.20909 

0.22491 

0.23418 

118.39 

0.25443 

75.910 

25 

0.10737 

0.10819 

0.11075 

117.64 

0.11599 

81.963 

0.20909 

0.21916 

0.22675 

118.07 

0.24232 

79.985 

50 

0.10737 

0.10830 

0.10965 

114.09 

0.11215 

90.156 

0.20909 

0.21598 

0.22022 

114.83 

0.22757 

88.475 

100 

0.10737 

0.10761 

0.10823 

111.36 

0.10949 

95.064 

0.20909 

0.21199 

0.21403 

111.92 

0.21772 

93.938 

200 

0.10737 

0.10735 

0.10769 

110.46 

0.10829 

97.569 

0.20909 

0.21011 

0.21124 

110.76 

0.21297 

96.925 


Table 5: Calculations relating to and R s (t) 
n = 20, c = 25, 6> = 0.8, t = 25, k = 2, m = 8 


r 

Col.l 

Col.2 

Col. 3 

Col.4 

Col. 5 

Col.6 

Col.7 

Col.8 

Col. 9 

Col. 10 

Col. 11 

Col. 12 

0 

0.00377 

0.00552 

0.00556 

101.87 

0.00471 

101.81 

0.00039 

0.00182 

0.00185 

94.865 

0.00162 

88.625 

5 

0.01152 

0.01420 

0.01436 

106.003 

0.01366 

86.929 

0.00355 

0.00885 

0.00914 

113.05 

0.00904 

79.946 

10 

0.03518 

0.03885 

0.03976 

112.85 

0.04073 

78.387 

0.03009 

0.04675 

0.04916 

120.39 

0.05214 

72.434 

15 

0.10737 

0.10779 

0.11066 

119.41 

0.11756 

78.399 

0.20909 

0.21937 

0.22761 

118.72 

0.24815 

76.659 

20 

0.32768 

0.31915 

0.32740 

138.77 

0.34411 

89.196 

0.79548 

0.75086 

0.75481 

117.19 

0.79071 

121.971 




















































Table 6: Calculations relating to R(t) and R s (t) 


n = 20, r = 15, c = 25, t = 25, k = 2, m = 8 . 


6i 

Col.l 

Col. 2 

Col.3 

Col.4 

Col. 5 

Col.6 

Col. 7 

Col. 8 

Col.9 

Col. 10 

Col. 11 

Col. 12 

0.4 

0.00010 

0.00032 

0.00032 

101.01 

0.00012 

337.78 

3.07 x 10“ 7 

1.49 x 10" 5 

1.50 x 10~ 5 

107.31 

4.39 x 10“ 6 

491.22 

0.5 

0.00097 

0.00183 

0.00184 

103.22 

0.00099 

204.17 

2.65 x 10 -5 

0.00028 

0.00029 

120.32 

0.00013 

230.03 

0.6 

0.00605 

0.00823 

0.00831 

104.67 

0.00597 

136.80 

0.00099 

0.00367 

0.00367 

111.03 

0.00236 

162.16 

0.7 

0.02824 

0.03160 

0.03213 

110.33 

0.02798 

102.49 

0.01995 

0.03348 

0.03482 

118.42 

0.02875 

112.80 

0.8 

0.10737 

0.10815 

0.11123 

118.94 

0.10703 

90.47 

0.20909 

0.22019 

0.22923 

119.02 

0.21801 

92.17 

0.9 

0.34867 

0.34157 

0.34866 

141.31 

0.34973 

97.45 

0.82891 

0.78709 

0.78698 

120.64 

0.79895 

105.99 

0.93 

0.48398 

0.47655 

0.48212 

165.12 

0.48399 

105.66 

0.95725 

0.93552 

0.92884 

155.35 

0.94084 

117.34 

0.96 

0.66483 

0.66788 

0.66198 

198.66 

0.66456 

118.19 

0.99731 

0.99469 

0.99086 

355.23 

0.99487 

120.60 

0.99 

0.90438 

0.93525 

0.90791 

173.29 

0.90454 

119.21 

0.99999 

0.99999 

0.99996 

390.01 

0.99999 

121.53 




Table 7: Unbiased Estimator 


n = 10, r = 15, c = 25, 6 = 0.8, k = 2, m = 8 


Combination \ t 

20 

25 

30 

35 

40 

+1-1+1-1+1-1+1-1+1-1 

334142.2 

30824.96 

8892.633 

13735.52 

-9741.06 


37026.8 

2746.153 

87121.23 

-23797.94 

4468.285 

+1+1000000-1-1 

-78456.49 

-31337.49 

37854.66 

843.4048 

-17.16015 

+100000000-1 

50794.91 

-36774.63 

14240.32 

1158.73 

-2478.892 

Reliability 

0.32768 

0.10737 

0.03518 

0.01152 

0.00378 


Table 8: Unbiased Estimator 


n = 10, r = 15, c = 25, t = 25, k = 2, m = 8 


Combination 4- | 0 —*■ 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

0.93 

0.96 

0.99 

+1 — 1 + 1 — 1 + 1 — 1 + 1 — 1 + 1 — 1 

277.34 

1297.87 

-7990.63 

-25294.98 

29730.69 

466268.7 

-626848.4 

1314555 

2910383 

+!+!+!+!+!—!—!—!—!—! 

128.73 

603.68 

2339.78 

2376.59 

-12487.09 

-378417.4 

-551749.3 

5398500 

-532512.9 

+1+1000000-1-1 

78.22 

-118.91 

3483.52 

21081.09 

-60072.55 

-107717.3 

-605752.3 

1757503 

-13004957 

+100000000-1 

132.79 

34.93 

-941.95 

3206.99 

-5265.09 

120432.4 

-695584.4 

530686.8 

-9375230 

Reliability 

0.000105 

0.000977 

0.006047 

0.028248 

0.107374 

0.348678 

0.483982 

0.664833 

0.904382 





Table 9: Calculation of 


n = 20, r = 15, c = 25, t 

6 Reliability Unbiased Estimator 




Mean 

Variance 

(LCL, UCL) 

0.4 

0.000104 

0.000123 

1.524 x 10- 7 

(0.000642, 0.000888) 

0.5 

0.000977 

0.001016 

4.575 x 10' 6 

(0.003176, 0.005208) 

0.6 

0.006046 

0.005932 

5.819 x 10” 5 

(0.009019, 0.020883) 

0.7 

0.028247 

0.028009 

0.000543 

(0.017681, 0.073699) 

0.8 

0.107374 

0.107834 

0.002879 

(0.002658, 0.213009) 

0.9 

0.348678 

0.349326 

0.006850 

(0.187109, 0.511543) 

0.93 

0.483837 

0.483837 

0.006413 

(0.326884, 0.640790) 

0.96 

0.664832 

0.665289 

0.004414 

(0.535067, 0.795511) 

0.99 

0.904382 

0.904440 

0.002129 

(0.813998, 0.994883) 


Interval and CP of R(t.) 


25, Repetation = 10000 

MLE 


CP 

Mean 

MSE 

LCL, UCL 

CP 

0.9698 

0.000309 

4.746 x 10“ 7 

(0.001041, 0.001659) 

0.9862 

0.9636 

0.001860 

8.889 x 10“ 6 

(0.003983, 0.007704) 

0.9830 

0.9450 

0.008182 

7.887 x 10“ 5 

(0.009225, 0.025589) 

0.9701 

0.9500 

0.031625 

0.000556 

(0.014591, 0.077841) 

0.9559 

0.9599 

0.108912 

0.002606 

(0.008854, 0.208971) 

0.9523 

0.9520 

0.341149 

0.006688 

(0.180857, 0.501442) 

0.9487 

0.9518 

0.476420 

0.006768 

(0.315177, 0.636989) 

0.9569 

0.9490 

0.668659 

0.005214 

(0.527134, 0.810185) 

0.9626 

0.8560 

0.935051 

0.002511 

(0.836834, 1.000000) 

0.9739 




Table 10: Efficiency of Rm with respect to R* M 

n = 10, r 2 = 5, 6>1 = 0.7, 6*2 = 0.8, R = 0.8543644 


Cl C 2 

15 

20 

25 

30 

10 

170.3204 

168.5565 

162.8126 

184.5733 

15 

123.3550 

119.4340 

121.8824 

124.4574 

20 

107.4592 

108.5157 

110.6393 

108.1562 

25 

106.2693 

104.6287 

105.2644 

103.2766 


Table 11: Efficiency of Rm with respect to R* M 

r, =5, r 2 = 10, 9 i = 0-8, 0 2 = 0.7, R = 0.8212655 


Cl c 2 

15 

20 

25 

30 

10 

147.774 

154.9791 

157.9302 

162.8920 

15 

114.8463 

118.9883 

118.0901 

121.116 

20 

110.571 

109.2363 

106.1009 

111.0485 

25 

103.5806 

104.4578 

100.1903 

103.9806 


Table 12: Efficiency of Rm with respect to R* M 

n = 10, r 2 = 5, gi = 0.7, 6*2 = 0.8, = 0.2234182 


Cl c 2 

10 

15 

20 

25 

15 

147.8276 

113.0031 

108.1288 

103.5886 

20 

159.4114 

115.1075 

105.5036 

104.5615 

25 

150.5815 

117.2885 

104.9115 

102.4807 

30 

140.7334 

117.5558 

107.1203 

102.5214 


Table 13: Efficiency of Rm with respect to R* M 

n = 10, r 2 = 5, 6*1 = 0.8, 6*2 = 0.7, 1? = 0.07639545 


Cl c 2 

10 

15 

20 

25 

15 

133.3005 

109.2876 

104.7854 

103.7044 

20 

140.5820 

108.7528 

101.8761 

100.2199 

25 

144.2547 

107.1663 

100.6334 

101.4030 

30 

146.9855 

114.6641 

102.8547 

101.3404 






































































































Table 14: MSEs of Estimator R 


r 2 |ri 

5 

10 

15 

20 


0.909091 

9.09091 xlO" 6 

9.09091 xlO" 11 

9.09091 xlO" 16 


0.920557 

0.000162 

1.37222xl0" 7 

1.55348xl0" 9 

5 

0.914516 

0.343349 

0.362465 

0.344545 


0.005747 

4.09776 xlO" 7 

8.34178xl0 -11 

3.25312xl0" 16 


0.006513 

0.280415 

0.300813 

0.284422 


0.999999 

0.909091 

9.09091 xlO" 6 

9.09091 xlO" 11 


0.9999468 

0.918347 

0.000160 

2.19983xl0" 7 

10 

0.999996 

0.91221 

0.350020 

0.347819 


4.47059 xlO" 08 

0.005862 

8.72553xl0" 7 

1.81692xl0" 12 


6.68094xl0" n 

0.006658 

0.288780 

0.283411 


1 

.999999 

.909091 

9.09091 xlO 6 


1 

0.999964 

0.916321 

0.000191 

15 

1 

0.999998 

0.909904 

0.331967 


2.63535xl0" 13 

2.19439xl0" 8 

0.005932 

4.68203 xlO" 7 


8.26455xl0" 23 

2.72526xl0" n 

0.006783 

0.271694 


1 

1 

0.999999 

0.909091 


1 

0.999999 

0.999964 

.914222 

20 

1 

1 

0.999998 

0.907743 


1.05866xl0" 18 

9.20165xl0" 13 

2.28706 xlO" 8 

0.006325 


1.77863xl0“ 32 

8.26455xl0" 23 

2.72526xl0" n 

0.007260 


26 






























Table 15: MSEs of Estimator R 


r 2 ri 

5 

10 

15 

20 


0.666667 

0.020833 

0.000651 

2.03450xl0 -5 


0.679177 

0.026838 

0.002103 

0.000206 

5 

0.671748 

0.021386 

0.001522 

0.001859 


0.012088 

0.000737 

2.33284xl0“ 5 

5.11005xl0" 7 


0.012973 

0.000982 

3.90682x 10 -5 

0.000701 


0.989583 

0.666667 

0.020833 

0.000651 


0.984516 

0.668551 

0.026998 

0.001960 

10 

0.989622 

0.660862 

0.021230 

0.001143 


0.000364 

0.011225 

0.000801 

1.61857xl0 -5 


0.002692 

0.012219 

0.000720 

2.13649xl0" 5 


0.999674 

0.989583 

0.666667 

0.020833 


0.998673 

0.984826 

0.669140 

0.0265948 

15 

0.999563 

0.990105 

0.661459 

0.020817 


1.50795xl0 -5 

0.000312 

0.012265 

0.000715 


6.954222 xl0~ 6 

0.000227 

0.013318 

0.000845 


0.999989 

0.999674 

0.989583 

0.666667 


0.999867 

0.998716 

0.984689 

0.670912 

20 

0.999988 

0.999662 

0.989835 

0.663109 


2.96125xl0" 7 

8.48119xl0" 7 

0.000345 

0.012157 


2.92646 xl0“ 8 

1.94865xl0 -6 

0.000254 

0.013208 


27 






























Table 16: MSEs of Estimator R 


r 2 |?’i 

5 

10 

15 

20 


0.238095 

7.61904xl0" 5 

2.43809xl0~ 8 

7.80190xl0 -12 


0.239828 

0.000384 

3.17275x 10“ 6 

7.44942xl0- 8 

5 

0.242251 

0.028820 

0.031522 

0.029636 


0.009815 

1.03719x 10 -6 

4.49468x 10 _1 ° 

6.71449 xlO" 13 


0.007733 

0.005539 

0.006393 

0.005711 


0.750339 

0.238095 

7.61905x 10~ 5 

2.43809x10-8 


0.753009 

0.236090 

0.003744 

3.40036x 10" 6 

10 

0.742105 

0.239565 

0.033643 

0.032150 


0.010630 

0.008706 

9.65263x 10” 7 

4.16401 xlO" 10 


0.012344 

0.006546 

0.006595 

0.006113 


0.918191 

0.750339 

0.238095 

7.61905x 10” 5 


0.918237 

0.756132 

0.240675 

0.000352 

15 

0.920495 

0.746061 

0.243698 

0.032421 


0.003296 

0.012589 

0.009720 

8.85057xl0- 7 


0.004096 

0.014265 

0.007625 

0.006236 


0.973193 

0.918191 

0.750339 

0.238095 


0.968532 

0.919645 

0.767043 

0.240174 

20 

0.973493 

0.922001 

0.757469 

0.242044 


0.000989 

0.003383 

0.011126 

0.009538 


0.001103 

0.004192 

0.012577 

0.007625 


28 






























Table 17: MSEs of Estimator R 


r 2 |?’i 

5 

10 

15 

20 


0.526316 

0.310784 

0.183515 

0.108364 


0.526810 

0.299460 

0.174113 

0.103238 

5 

0.526508 

0.308696 

0.182257 

0.105726 


0.017799 

0.012224 

0.007922 

0.0043331 


0.015585 

0.011825 

0.008655 

0.005117 


0.720294 

0.526316 

0.310784 

0.183515 


0.726634 

0.529795 

0.305464 

0.175418 

10 

0.716990 

0.529132 

0.314470 

0.183407 


0.011833 

0.016580 

0.013347 

0.008355 


0.011819 

0.014537 

0.012886 

0.009171 


0.834836 

0.720294 

0.526316 

0.313784 


0.840993 

0.722299 

0.530707 

0.305661 

15 

0.833705 

0.712836 

0.529979 

0.314475 


0.007389 

0.012591 

0.015473 

0.013614 


0.008308 

0.012569 

0.013530 

0.013235 


0.902473 

0.834836 

0.720294 

0.526316 


0.903279 

0.840593 

0.726930 

0.525431 

20 

0.900223 

0.833145 

0.717270 

0.525064 


0.003855 

0.006583 

0.012184 

0.017255 


0.004660 

0.007377 

0.012206 

0.015104 


29 






























Table 18: MSEs of Estimator R 


r 2 |?’i 

5 

10 

15 

20 


0.952381 

0.321076 

0.102261 

0.033509 


0.954791 

0.300429 

0.107198 

0.038337 

5 

0.953994 

0.313513 

0.107466 

0.033587 


0.001314 

0.012567 

0.005071 

0.001269 


0.001318 

0.013776 

0.006272 

0.001443 


0.999985 

0.952381 

0.312076 

0.102261 


0.999870 

0.952055 

0.296181 

0.103630 

10 

0.999989 

0.951173 

0.309150 

0.103415 


1.49867xl0 -7 

0.001352 

0.012463 

0.004663 


6.16543xl0" 9 

0.001373 

0.013512 

0.005783 


1 

.999985 

.952381 

.312076 


0.999998 

0.999865 

0.952267 

0.302442 

15 

1 

0.999983 

0.951499 

0.315536 


1.01544xl0 _1 ° 

2.91485xl0 -7 

0.001478 

0.012921 


2.37772xl0 -17 

4.87421 xlO" 8 

0.001480 

0.014082 


1 

1 

0.999985 

0.952381 


1 

.999998 

0.999858 

0.952217 

20 

1 

1 

0.999986 

0.951531 


4.40539xl0 -13 

1.37763xl0 _1 ° 

1.86115x 10 -7 

0.001436 


2.43492xl0 -24 

2.37772 xlO" 17 

8.33376x 10 -9 

0.001436 
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Table 19: Calculation of Confidence Interval and CP of R 


n i = n 2 = 10, 9 1 = $2 = 0.9, Repetation = 10000 


(ri, r 2 ) 

Reliability 

Unbiased Estimator 

MLE 

Mean 

Variance 

(LCL, UCL) 

CP 

Mean 

MSE 

LCL, UCL 

CP 

(5, 20) 

0.108363 

0.108382 

0.005508 

(0.037085, 0.253851) 

0.9540 

0.104635 

0.004648 

(0.028995, 0.238266) 

0.9372 

(5, 15) 

0.183515 

0.184696 

0.009011 

(0.001353, 0.370744) 

0.9650 

0.176514 

0.008193 

(0.000888, 0.353917) 

0.9510 

(5, 10) 

0.310784 

0.309366 

0.012371 

(0.091370, 0.527362) 

0.9540 

0.300321 

0.012822 

(0.078386, 0.522256) 

0.9592 

(5, 5) 

0.526316 

0.524415 

0.014289 

(0.290129, 0.758701) 

0.9504 

0.524652 

0.016296 

(0.274403, 0.774854) 

0.9627 

(10, 5) 

0.720294 

0.721299 

0.011505 

(0.511075, 0.931525) 

0.9545 

0.730976 

0.011572 

(0.520138, 0.941814) 

0.9560 

(15, 5) 

0.834836 

0.832709 

0.008217 

(0.672409, 1.000000) 

0.9629 

0.840145 

0.007324 

(0.672409, 1.000000) 

0.9470 

(20, 5) 

0.902472 

0.902443 

0.004755 

(0.767295, 1.000000) 

0.9545 

0.905403 

0.003967 

(0.781953, 1.028853) 

0.9545 













